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Creation of scalar and Dirac particles in the presence of a time varying electric field in 

an anisotropic Bianchi I universe 
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In this article we compute the density of scalar and Dirac particles created by a cosmological 
anisotropic Bianchi type I universe in the presence of a time varying electric field. We show that 
the particle distribution becomes thermal when one neglects the electric interaction. 


I. INTRODUCTION 

During the last decades a great deal of effort has been made in understanding quantum processes in strong fields. 
Quantum field theory in the presence of strong fields is in general a theory associated with unstable vacua. The 
vacuum instability leads to many interesting features; among them particle creation is perhaps the most interesting 
nonperturbative phenomenon. 

In order to analyze the mechanism of particle creation in cosmological backgrounds we have at our disposal different 
techniques such as the adiabatic approach |f|,^| , the Feynman path integral method ||, the Hamiltonian diagonalization 
technique [||j|], as well as the Green function approach 0 

The study of quantum processes in the presence of electromagnetic fields in curved backgrounds has almost been 
restricted to constant electric fields. The density of particles created by an intense electric field was first calculated 
by Schwinger |Q and different authors have discussed this phenomenon in curved spaces. The technical difficulties 
associated with solving the Klein-Gordon or Dirac equations in time dependent electromagnetic fields have reduced 
the number of configurations studied to a few simple solvable models ji|. A further difficulty arises when one deals 
with cosmological backgrounds with singularities [g|| . Here one cannot apply the standard adiabatic approach and 
other equivalent methods developed in the literature. It is worth mentioning that the presence of primordial electric 
fields could enhance the particle creation mechanism and also produce deviations from the thermal spectrum. 

In order to discuss particle creation in an anisotropic universe, we consider the cosmological scenario associated 
with the metric 

ds 2 = — dt 2 + t 2 (dx 2 + dy 2 ) + dz 2 (1) 

The line element (]!]) presents a singularity at t = 0, The scalar curvature is R = 2 /t 2 , and consequently, the adiabatic 
approach Jj| cannot be applied in order to identify particle states. Using the Hamiltonian diagonalization method, 
Bukhbinder and coworkers 100 were able to compute the density of particles created by the background associated 
with the metric (0), and also in Ref. ||T^] this result was extended to include a time dependent electric field. 

The purpose of the present article is to discuss the production of scalar and Dirac particles in the anisotropic 
cosmological background associated with the line element (|lj) in the presence of a time dependent homogeneous 
electric field. In order to compute the rate of particle creation we apply a quasiclassical approach that has been used 
successfully in different scenarios 00- The idea behind the method is the following. First, we solve the relativistic 
Hamilton-Jacobi equation and, looking at its solutions, we identify positive and negative frequency modes. Second, 
we solve the Klein-Gordon and Dirac equations and, after comparing with the results obtained for the quasiclassical 
limit, we identify the positive and negative frequency states. The paper is structured as follows. In Sec. II we solve 
the relativistic Hamilton-Jacobi equation and compute the quasiclassical energy modes. In Sec. Ill we solve the 
Klein-Gordon and Dirac equations and obtain the density of scalar and Dirac particles created. The discussion of the 
results and final remarks are presented in Sec. IV. 
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II. COMPUTATION OF THE QUASICLASSICAL SOLUTIONS 


The relativistic Hamilton-Jacobi equation can be written as 

s “'(S;- eA ‘)(Si- eAf ) +m2=0 (2) 

where here and elsewhere we adopt the convention c = 1 and h = 1. 

The vector potential A M 


& = j6% (3) 

corresponds to a time dependent homogeneous electric field Ck/t 2 . The invariants 

If^f iu/ = ~, f^f; i/ = o (4) 

indicate that there is no magnetic field. Since the line element (|l]) and the vector potential (|3]) do not depend on the 
space variables we can write the function S as follows: 


S = k ■ r + F(t). 


(5) 


Substituting Eq. (j^) into the Hamilton-Jacobi equation one obtains 


k 2 x + k 2 y 
t 2 


/T eG.2 , dF , n 2 r\ 

( ‘*-T) - ( a-> +m =°- 


The solution of (|(]) has the following asymptotic behavior 

lim F(t) = ±^k 2 x + k 2 y +e 2 C 2 log* 


( 6 ) 

(7) 


$ = e iS _> C'ot ± V fc » +fc ? +e2c2 


( 8 ) 


as t 


0, and 


_ UpC 

lim F(t) = ±y/k 2 + m 2 t T , 1 = log t , (9) 

t^+oo _|_ J[,f2 

T -- ■■ 

$ = e iS —> C'ie ±i V /fc ^ +m2t t \f k l+ M2 (10) 

as t —> oo. The positive or negative frequency mode is selected depending on the sign of the operator idt . Positive 
and negative frequency modes will have, respectively, positive and negative eigenvalues; thus in Eqs. ( 0 )-© the 
upper signs are associated with negative frequency modes and the lower signs correspond to positive frequency states. 
After making this identification we proceed to analyze the solutions of the Klein-Gordon and Dirac equations in the 
cosmological background (Q) . 


III. PARTICLE PRODUCTION 
A. Scalar particles 

The covariant Klein-Gordon equation has the form 

g af3 (V Q - ieA a ) (V^ - ieA 0 ) $ - (m 2 + £R)3> = 0 (11) 
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where V Q is the covariant derivative, R is the scalar curvature and £ is a dimensionless coupling constant which takes 
the value £ = 0 in the minimal coupling case and £ = 1/6 when a conformal coupling is considered. The scalar wave 
function $ is normalized according to the inner product fl^] 


($i,$2>= f ($20^ - $^$ 2 )V ZI 9ds 11 

J a 


( 12 ) 


where a is an arbitrary spacelike hypersurface. Since the metric (Jlj) is stationary we choose a hypersurface orthogonal 
to the timelike vector = <5 q . Since Eq. © commutes with each of the components of the linear momentum 
p = — zV, the wave function <f> can be written as 


$ — i _1 A(t)e i ( fc * x + k yy+ k z z )' 


(13) 


Substituting Eq. (^3|) into Eq. © we reduce the problem of solving Eq. © to that of finding a solution of the 
following second order differential equation: 


cP_ 

dt 2 


k 2 + k 2 + e 2 C 2 + 2£ „ 2 k z eC 2 


+ k, — 


t 2 ' * t 

The solution can be expressed in terms of Whittaker functions fis| , |l9| 

A = £ 1 M M (p) + £ 2 %(p) 


+ m 2 A (t) = 0 


(14) 


(15) 


where C 2 and C 2 are arbitrary constants, and p, k and p are given by the expressions 


_ / ] av. P n 

p = -2 i^/k 2 z +m 2 t, p = +kl + k 2 + e 2 C 2 +2 £,k = - ^ = (16) 

V Q -v / rv z 7TI 

Looking at the asymptotic behavior of M k ^{z) as z —> 0 

M k ^z)~e-*l 2 z x l 2 +», (17) 

we find that, according to (j|), and using the fact that all the coefficients in Eq. (0) are real, the positive and negative 
frequency solutions (|l5[) at t = 0 are 

A+ = Co + M k , fi {p) i Aq = (Co + M k ^(p))* = Co + (-l)" /i+1/2 M fc ,_ /2 ( /9 ) (18) 

where C+ is a normalization constant. Analogously, looking at the behavior of W k . /t (z) as \z\ —* oo 

W k Jz) - e~ z ' 2 z k , (19) 


we see that the corresponding positive and negative frequency modes as \p\ —> oo are 


A+ = £+W fclA (p), A" = C" W- k , fi {-p) 


( 20 ) 


where C+, and C^ are normalization constants. The positive frequency mode at \p\ —> oo can be expressed in terms 
of the positive Aj and negative A(^ frequency modes via the Bogoliubov transformation 

A+ =aA++/3A 0 ”. (21) 


Since the Whittaker function W k ^(z) can be expressed in terms of M k ^(z) as follows: 


W k>lt (z) = 2m) , M kilt (z) 


T(2p) 


Y{\- p-k) 

we obtain that a and /3 are given by the expressions 

C r(-2 p) 


T(4 + M - k ) 


Affc,—/x(^), 


a = 


<d" r(s - v - k) 


, (3 =-i 


r(2/2) 

€t r H + p-k) 


exp( —7T \p\) 


( 22 ) 


(23) 


then 
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(24) 


p 

2 

r(i-/x-fc) 

a 


r(| + ji — k) 


exp(—2?r |/Lt|) 


The computation of the density of particles created is straightforward from Eq. (^4j) and the normalization condition 
of the wave function which, for scalar particles, means that the Bogoliubov coefficients satisfy the relation 


\a\ 2 -\8\ 2 = l. 


Taking into account that 


r (^ +w) 


cosh ny ’ 


we find that the density of particles created is 

2 exp(7r(|/2| + |fc|))cosh(7r(|/i| + |fc|)) exp(-27r |/z|) 

71=1/31 = -sinh(2^|)-' 


(25) 


(26) 


(27) 


It is worth mentioning that, thanks to the relations (^) and (|24|), we did not have to compute the normalization 
constants (£+, C T, a nd £+. When the electric field is switched off, the parameter k is zero and the density of 
particles created (p7|) reduces to a thermal distribution ||: 


n = 


exp 27 T\ k\ + 2£ - i - 1 


(28) 


B. Dirac particles 

In order to compute the density of spin 1/2 particles created, we proceed to solve the Dirac equation in the 
cosmological background (|l!) in the presence of the electric field (j|). 

The covariant Dirac equation in curved space in the presence of electromagnetic fields can be written as follows 


- ieAn) + m}T = 0 


(29) 


where the curved gamma matrices satisfy the anticommutation relation {y^y 1 '} , = 2 g( tv and the spinor connec¬ 


tions T m are 


where 


^ g\o 


( db 0 \ 

I v I n a _ r“ 
\dx») a e ^ 




s = ^{1 ~l v l )■ 


The matrices b 0 , a^ establish the connection between the curved y^ and Minkowski y^ Dirac matrices as 

7 /. = b“%, 7 m = a$7 0 

and 

tV+tY = ‘2r] Xv . 

Since the line element ( li) is diagonal, we choose to work in the diagonal tetrad 

7 m = \/|5^ Al |7 A \ no sum. 

Substituting Eq. (]34|) into Eq. (pOl) we obtain 


(30) 

(31) 

(32) 

(33) 

(34) 
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( 35 ) 


Ti = \fi\ r 2 = r 3 = o, r 0 = o, 


and the Dirac equation (|29|) takes the form 


~o d l,.i 9 ~2 d \ -3/ & ieC , 

1 ai + l {l di + 1 ai ) + 1 , + ”d*° = tl 


(36) 


where To = tT. The factor t was introduced in order to cancel the contribution due to the spinor connections 
Equation (pfl) can be written as a sum of two first order commuting differential operators as follows [pO|| 


(A'-i + AT 2 )T = 0, 


(37) 


k 2 § = k<s> = —k\$, 

where the spinor $ is related to To via the expression 

7 3 7°T 0 = T, 

and k is a separation constant. The operators K\ and iv 2 read 


K, = t 




(38) 


(39) 


(40) 


dx 


dy 


^=(7 1 i-+7^)7 3 7°. 


(41) 


Since Eq. (|3(3|) commutes with —zV, the spinor $ can be written as $ = $o exp(z(fc x :r + k v y + k z z )). Choosing to 
work in the following representation of the Dirac matrices: 


7° = 


—ZCT 1 0 

0 ZCT 1 


1,7 = 


0 i 

-i 0 


7 2 = 


a 2 0 

0 -a 2 


2,7 = 


CT 3 0 

0 -cr 3 1 ’ 


(42) 


we see that the equation k 2 T = M> helps us to determine the relation between the components of the bispinor <f> 


To = 


$i 

<j> 2 


I 

\ ik v —K 


where 


k = i^Jk 2 + ky. 

Using the representation (|||) we find that Eq. (|o|) reduces to the system of equations 


cr 3 -^- + o' 2 (k z - -j-) + me 2 + 7 $i = 0 , 


(43) 


(44) 


(45) 


j-er 3 -^ + cj 1 {k z - ^-) + ma 2 + ^ j T 2 = 0. (46) 

Taking into account the structure of the bispinor (^|), it is not difficult to see that Eq. ([It]) is equivalent to Eq. ( |45| ) 
and consequently we have reduced the problem of solving Eq. (|||) to that of finding the solution of Eq. ([l5|) . 

In order to solve Eq. ( |45| ) we proceed as follows. Since the structure of Eq. ( ffo| ) resembles the one obtained after 
solving the hydrogen atom in spherical coordinates |f22| , [23|| and we are interested in identifying positive and negative 
frequency modes, it would be desirable to have solutions of Eq. (fl5|), as in the scalar case, in terms of single Whittaker 
functions. For this purpose, we apply to the spinor $i the linear transformation T 
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T$ i = F = 


with 

The components of the spinor F 

dfi , 

- /1 

dp 


T = 


\J im. + k z — \/ im — k z 


i /im + k z \J—im — k z 
[|) satisfy the system of equations 


-Cek z 


- 


m 


+ ?\-h 


IeCm 


2 2 


\f I W^ 


- k] = 0, 




(47) 


(48) 


(49) 


dh f I _ 

^ dp i y/— k z — m 2 2 


- Cefe + S)+A 


IeCm 


y/—k& — m 2 


fc = 0, 


(50) 


whose solutions can be expressed in terms of the Whittaker functions Wfc, M (;z) and Mk^{z) fl8| , |l9|] . 

In order to construct the positive and negative energy frequency modes, we use the asymptotic behavior of the 
hypergeometric functions © and © and compare the solutions of Eqs. (^) and (|50]) with Eqs. (||) and ©), 
obtained after solving the Hamilton-Jacobi equation. Using this procedure, we identify the positive frequency mode 
as t —> Too as 


F + = 

r -(-OO 


( £oo + P 1/2 W 1[ /,. 2 

2 yj-kl-m 2 

, } Ce ,n - P~ l/2 Wl _ eC k ■■ (p) 

\ \J~ k z-rn - 2 -y/-fc 2 -m 2 


with 


p = V fc 2 — e 2 C 2 = * y + e 2 C 2 = ip 

£+ is a normalization constant according to the inner product 

(4>i,4> 2 )= / $ 1 7 A1$ 2 V-gda^ 

J a 

where a is an arbitrary spacelike hypersurface. Analogously, we have that as t —> 0 the positive mode is 


(51) 


(52) 


(53) 


/ 


F + — 
r o — 


£ 0 +p- 1 / 2 M i 


P 1 + eCt, _ „(p) 


yj-k 2 -m 2 ' 




\ 


f+M 


V 


» 


(54) 




where £+ is also a normalization constant. Using the asymptotic behavior of the Whittaker function Mk^{z) ©) 
and the normalization condition (E3), we obtain the result for the negative frequency mode F 0 as t —> 0: 


Fo = 


( 1/2 M l , eCkz (p) 

ecfcz -/X 

e^g+p- 1 / 2 !^,, ec, „(p) 


\ 


V 


(55) 


/ 


V- fe i— 2 ’ 

The density of particles created can be calculated with the help of the Bogoliubov transformation 

F+=aJ’ 0 ++^F 0 - 


(56) 


+ r (~ 2 u0 F + 

+oc T(i-ip-e) 0 


T(2ip) 


r(| + v-e) 


0 -7T/^ 


F f 


o > 


(57) 
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where € = —i H- = = = . Using the property of the Gamma function 

z y —kj. —m 2 


l r (*y)l 2 = 


y sinh Try 


we obtain 


eCk z 


sinh7r(/i t , 
0 -2-Kfi _ V k *+r 


eCk z 
\Jk\-\-v 


sinh7r(/i — 


eCk z 

\A T+~ r - 


eCk z 


\J kl +n 


From Eq. 


and taking into account the normalization condition for Dirac particles 


= 1 


we find that the density of Dirac particles created is 


exp(—27r/i) ( fi — 


eCk z 

\/k 2 +n 


sinh 7 r(/i ■ 


eCk z 

y/kz+n 


exp(—27 tu) (u - , cCjt ' ; ) sinh7r(/i H — eCk = ) + (h H— , eCfr - ) sinh7r(/i- f ck * ) 

y yk 2 z +m 2 J ykl+m 2 ' y yk 2 +m 2 J ^ yk 2 z +m 2 ’ 

In the absence of an electric field, the density of particles created © becomes thermal: 

1 

TL == -. 

exp27r |fcj_| + 1 


(58) 


(59) 


(60) 


(61) 


(62) 


IV. DISCUSSION OF THE RESULTS 


From the expressions (27) and (61) we observe that, in the presence of the electric field (y), the distribution of 
particles created is not thermal and it strongly depends on the strength eC. For very strong electric fields we have 
that Eq. (E?j) reduces to 


n ~ eM-^\Jkl+e>C* + 2d-\ + -^==1 (63) 

an analogous result can be obtained for spin 1/2 particles. Expression ( |63| ) corresponds to a thermal distribution 
with an effective mass and a chemical potential proportional to eC . Equation ( |63|) shows that the presence of 
strong electric fields contributes significantly to the creation of particles. The results obtained in this article show 
that the quasiclassical approach permits one to compute positive and negative frequency modes even when spacetime 
and electromagnetic sources are not static, and encourage us to study quantum effects in more realistic anisotropic 
scenarios. 
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